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Abstract 

The influence of quantum degeneracy on the performance of a Stirling refrigeration cycle is 
investigated, based on the equation of state of an ideal Fermi gas. The inherent regenerative 
losses and the coefficient of performance (COP) of the cycle are calculated. It is found that, 
under the condition of strong gas degeneracy, the COP of the cycle in the first approximation 
is a function only of the temperatures of the heat reservoirs, while under other conditions, the 
COPs of the cycle depend on the temperatures of the heat reservoirs and other parameters of 
the cycle. The results obtained here reveal the general performance characteristics of a Stirling 
refrigeration cycle having a Fermi gas as its working substance. 

© 2002 Elsevier Science Ltd. All rights reserved. 
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1. Introduction 

Recently, quantum-mechanical cycles have become interesting research subjects. 
The performances of quantum heat engines and refrigerators have been discussed 
[1-7]. Like classical cycles, quantum cycles may have different working substances. 
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Nomenclature 

Cy Heat capacity at constant volume (J) 

F(z) Correction factor 

f(z) Fermi function 

g Weight factor 

h Planck’s constant (Js) 

k Boltzmann’s constant (J/K) 

m Rest mass of an ideal Fermi particle (kg) 

N Total number of particles of an ideal Fermi gas 

n Number density of of an ideal Fermi gas (m -3 ) 

P Pressure (Pa) 

Qbc Amount of heat exchange between the working substance and the 
regenerator for a constant F L process (J) 

Qda Amount of heat exchange between the working substance and the 
regenerator for a constant Vh process (J) 

Qw Amount of heat exchange between the working substance and the hot 
reservoir for a the high temperature process (J) 

Q L Amount of heat exchange between the working substance and the cold 
reservoir for a low temperature process (J) 

Q 1 Amount of heat exchange for a isochoric process (J) 

r v Volume ratio of two isochoric processes 

S Entropy of the gas (J/K) 

T Temperature of the gas (K) 

T f Fermi temperature (K) 

T\\ Temperature of the hot reservoir (K) 

T\ Temperature of the cold reservoir (K) 

U Internal energy (J) 

V Volume (in 3 ) 

V H Maximum volume of the gas (m 3 ) 

Vp Minimum volume of the gas (m 3 ) 

v Average volume that an ideal Fermi particle occupies (m 3 ) 

vh Maximum average volume that an ideal Fermi particle occupies (m 3 ) 

v>L Minimum average volume that an ideal Fermi particle occupies (m 3 ) 

W Work input per cycle (J) 

z Fugacity of the gas 

A Q Regenerative loss (J) 

s Coefficient of performance 

e c Coefficient of performance of a Carnot refrigeration-cycle 

X Mean thermal wavelength (m) 

// Chemical potential of the gas (J) 
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For example, the working substance in quantum cycles may be harmonic oscillator 
systems, spin systems, ideal quantum gases, and so on. Thus, it is necessary to 
study the influence of a specific working substance on the performance of a quan¬ 
tum cycle. 

According to the theory of classical thermodynamics, a Stirling refrigeration cycle 
working with an ideal classical gas [8,9], a homogeneous ferromagnetic material 
[10,11], or a homogeneous ferroelectric material [12] may possess the condition of 
perfect regeneration through the use of a reversible regenerator. Its coefficient of 
performance (COP) is equal to that of a reversible Carnot refrigeration cycle. 
Flowever, when the working substance is an ideal quantum gas, this conclusion is 
not valid because the quantum degeneracy of the gas particles becomes important 
[13-15]. Therefore, it is of significance to analyze how the performance of the Stir¬ 
ling refrigeration cycle working with an ideal quantum gas is affected by the quan¬ 
tum degeneracy of the working substances. 

In the present paper, based on the equation of state of an ideal Fermi gas, the 
inherent regenerative losses and the COP of a Stirling regeneration cycle using an 
ideal Fermi gas as the working substance are analyzed. Under the conditions of 
quantum degeneracy, the ratio of the COP of the cycle to that of a reversible Carnot 
refrigeration cycle versus the temperature ratio of two heat reservoirs and the 
volume ratio of two isochoric processes are examined. 


2. Thermodynamic properties of an ideal Fermi gas 


Based on quantum statistics, the expressions of the pressure and the number den¬ 
sity of an ideal Fermi gas are given by [13] 


P = 


gkT f < \ 
A 


(i) 


and 


N 1 8 f m 

n=y=- = gfip.V ) ’ ( 2 ) 

respectively, where g is a weight factor that arises from the “internal structure” of 
the particles (such as spin), A = h/(2nmkT) l/ 2 is the mean thermal wavelength of the 
particles, h is Planck’s constant, m is the rest mass of a particle, k is Boltzmann’s 
constant, T is the gas temperature, N is the total number of the particles, V is the 
volume of the system, v is the average volume that a Fermi particle occupies, 
2 = exp (/j./kT) is the fugacity of the gas, fi is the chemical potential of the gas, 
1 00 x u ~ 1 d x 

f n (z) = - —:- : — is called the Fermi function, and F(«) is the Gamma func- 

F(«)Jo 2 _1 e v +l 

tion. By using Eqs. (1) and (2), the correction equation of state for an ideal Fermi 
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gas may be written as 
P = nkTF(z), 


(3) 


where F(z) is called the correction factor and it is defined as 


F(z) = 


hr 2(-) 

h/ih) 


(4) 


The internal energy and the entropy of an ideal Fermi gas are, respectively, given by 
[13] 

U=\ = 3 Nk TF(z) (5) 

2 fy 2 (z) 2 

and 


S — Nk 

zi 

= Nk 

1 

*| 

1 

5“ 

_ 1 


.2/3/2 (z) 


|2 J 


From Eq. (5), one obtains the heat capacity at constant volume as 

c ^(w)4 m h TF ^- 


(7) 


where F(T, v) is a function of temperature T and volume v. 

Using Eqs. (6) and (7), one can calculate the amounts of heat exchange in an iso¬ 
thermal and isochoric process. They are given by 


gif = T(S f - SO = NkT 


, [F(z f ) - F(z t )\ - (In z f - In z ; ) 


and 


( 8 ) 


Qit = 


C V (T, v)d T=-Nk[T f F(T { , v) - FF(T U v)], 


(9) 


where the indices i and /refer to the initial and final states, respectively. 

With the help of Eqs. (8) and (9), we can analyze the inherent regenerative losses 
and the COP of a Stirling refrigeration cycle working with an ideal Fermi gas. 


3. A Fermi-Stirling refrigeration cycle 

A Stirling refrigeration cycle using an ideal Fermi gas as the working substance is 
called a Fermi Stirling refrigeration cycle, which is composed of two isothermal and 
two isochoric processes. The entropy-temperature diagram of the cycle is shown in 
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Fig. 1. The entropy-temperature diagram of a Fermi Stirling refrigeration cycle. 


Fig. 1, where Q L and Q\\ are the amounts of heat exchange between the working 
substance and the heat reservoirs at temperature F L and r H during the two 
isothermal processes, Q bc and Q da are the amounts of heat exchange between the 
working substance and the regenerator during the two isochoric processes, and Vh 
and Fl are the maximum and minimum volumes of the gas system. All heats Q L , 
Q h , Q bc and g da are positive. In order to improve the performance of the cycle, a 
regenerator is often used in two isochoric processes. 

By using Eqs. (8) and (9), the amounts of heat exchange in the various processes 
of the Stirling refrigeration cycle may be expressed as 


Q l = 5 -NkT L [F(z A ) - F(z q )\ - NkT L (\n z d - In z c ), 

(10) 

6h = *NkT H [F(z a ) - F(z b )] - NkT H (\n z a - In z b ), 

(11) 

Qda = Cv(T, v H )dr = ^Nk[T H F(z a ) - T L F(z d )], 
id l 

(12) 

and 


[ h 3 

gbc = C V (T, v L )dF= -Nk[T H F(z b ) - T h F(z c )\, 

C ^ 

(13) 


respectively, where z a , z b , r c , and z d are the fugacity of the gas in the a, b, c, and d 
states of Fig. 1. 
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Fig. 2. The function TF{T , v) versus the temperature curves for two different volumes. The solid curve is 
obtained for a constant minimum average volume, vl = 10 -29 m 3 , while the dashed curve is obtained for a 
constant maximum average volume, i>h = 10 -28 m 3 . 


From Eqs. (10) (13), we obtain the work input per cycle as 
w = Qh — Ql + 2bc — (?da 

= NkT H [F(z. d ) - F(z b ) + lnr b - lnr a ] + NkT L [F(z c ) - F(z d ) + Inz d - lnz c ] 


(14) 

The curves of the function TF( T. v) varying with the temperature for two different 
volumes are shown in Fig. 2. Using Eq. (7) and comparing the slopes of the curves in 
Fig. 2, we can conclude that 


C V (T, vh) > C V (F, vl) 


(15) 


and 


A Q — Qbc ~ Qda 

= \Nk{T n [F(z h ) - F(z d )\ + F h [F(z d ) - F(z c )]} < 0, (16) 

where v H and v L are the maximum and minimum average volumes that an ideal 
Fermi particle occupies. Eq. (16) shows that the amount of heat exchange Q hc 
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transferred into the regenerator in one isochoric process is smaller than that of heat 
exchange Q da transferred from the regenerator in the other isochoric process. The 
inadequate heat in the regenerator per cycle can only be compensated from the hot 
reservoir in a timely manner, while refrigeration heat Q L is unvarying. If not, the 
temperature of the regenerator would be changed such that the regenerator would 
not be operated normally. From Eqs. (10) and (14), the COP of the cycle can be 
expressed as 


e = 


Ql 

w 


^ T L [F(z d ) - F(z c )] - 71(1 nz d - lnz c ) 

T H [F(z a ) - 7 7 (z b ) + lnz b - lnz a ] + T L [F(z c ) - F(z d ) + lnz d - lnz c ]' 


(17) 


By using Eqs. (16) and (17), the non-perfect regeneration characteristics and the 
COP of the cycle may be analyzed in detail. 


4. Three special cases 


(1) Under the condition of strong gas degeneracy and in the first approximation, 
we can derive the fugacity and the correction factor (a detailed derivation is given in 
Appendix A) as 


lnz = 


Ty(v) 

T 


ir T 

nfTiv) 


( 18 ) 


and 


F(T, v) = 


27V(v) 
5 T 


it 2 T 
67>(v) ’ 


(19) 


where 71(v) = 7)/v 2/3 is the Fermi temperature and D = (3/i 3 /7r)“ /3 /(8 km). 

By substituting Eqs. (18) and (19) into Eqs. (16) and (17), the inherent regen¬ 
erative loss and the COP of the cycle in the first approximation can be expressed as 


= jD Nk ( T H - r L 2 )('l 2/3 - V(T) < 0 


( 20 ) 


and 


T l 2T l 
T h ~T l T h + T l 


( 21 ) 


respectively, where e c = 71/(71 — 71) is the COP of a reversible Carnot refrigeration 
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cycle. In such a case, the COP is only a function of temperature and independent of 
any other parameters. 

(2) Under the condition of weak gas degeneracy, to a first approximation, we can 
derive the fugacity and the correction factor (a detailed derivation is given in 
Appendix B) as 

z = 4V2E/(T 3/2 v)[l + 2£'/(7’ 3/2 v)] (22) 


and 


F(T, v) = 1 +E 


1 

tw 


(23) 


where E = h 3 /2>2{n ink) 312 . 

By substituting Eqs. (22) and (23) into Eqs. (16) and (17), the inherent regen¬ 
erative loss and the COP of the cycle can be determined by 

A Q = \NkE(vp - Vh‘)(i/E jf - 1 /Tl' 2 ) < 0 (24) 


and 


r L ln(v H /v L ) - (1/2)(e/^ /2 )(1/v l - 1/vh) 


(T h - r L )ln(v H /v L ) + E(l/v L - 1 /v h )[(i/7h / 2 ) - (l/^ /2 )‘ ' 


(25) 


When 3 ELe gas is chosen as the ideal Fermi gas and Eq. (25) is used, one can generate 
the curves of the ratio of the COPs of the cycles, e/e c , varying with the temperature 
ratio, r=r H /CL, of two heat reservoirs for three different volume ratios, 
r v — vh/vl, of two isochoric processes, as shown in Fig. 3, where the minimum 
average volume that a particle occupies and the highest temperature are chosen to be 
vl = 10 _29 m 3 and 7 h = 177 K, respectively. For a constant maximum average 
volume value, vh = 10 _28 m 3 , or a constant minimum average volume value, 
vl = 10 _29 m 3 , Fig. 4 gives the curves of the ratio of the COPs of the cycles varying 
with the volume ratio of two isochoric processes, where the highest and lowest 
temperatures are given as 7 h = 177 K and 7) = 7)( i'i ) = 16.6 K, respectively. For 
a constant vl, the coefficient of performance increases with an increase of r v , while 
for a constant vh, the coefficient of performance decreases with an increase of r v . 
This implies that the COP of the cycle may be improved by increasing the maximum 
volume of the cycle. 

When the temperature is high enough, (1/2)^£ , /7 , ] 3/2 ^(1 /vl — l/vH)/ln(r v ) << 1 
and E( l/v L — l/v H ) (1 /7’ f l [ " 2 j — ^1 /T [ l/2 ^ /ln(r v ) << 1. Eq. (25) is simplified as 


T h ~T l 


Be¬ 


rn 
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Fig. 3. The ratio of the COP of a Fermi Stirling refrigeration cycle to that of a reversible Carnot refrig¬ 
eration cycle versus the temperature ratio curves for three different volume ratios of two isochoric pro¬ 
cesses. The curves are presented for 3 FIe. The minimum average volume and the highest temperature are 
given as vl = 10 _29 m 3 and Th = 177K, respectively. 


In this case, an ideal Fermi gas becomes an ideal classical gas and the working gas is 
always in the classical gas state throughout the cycle. The Stirling refrigeration cycle 
utilizing an ideal classical working gas may possess the condition of perfect regen¬ 
eration, so that its COP is equal to that of a reversible Carnot refrigeration cycle. 

(3) Under the condition of Th >> 7f(v l ) and C L << T F (v H ), the fugacity and 
correction factor may be expressed as 


z(T H , v) = AV2 E/(tK 2 v) [l + 2 E/(t% 2 v) 


(27) 


F(T H ,v) = l+E/(T% 2 v y 


lnz(T L , v) = 


T f (v) 

T l 


TT 2 T h 
12 7>(v)’ 


and 


F(T l , v ) = 


2C f (v) 
5 C L 


77 2 T L 

67p(v)' 


(28) 

(29) 


(30) 


respectively. By substituting Eqs. (27) (30) into Eqs. (16) and (17), the inherent 
regenerative loss and the COP of the cycle are determined by 
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Fig. 4. The ratio of the COP of a Fermi Stirling refrigeration cycle to that of a reversible Carnot refrig¬ 
eration cycle versus the volume ratio curves for 3 He. The solid curve is obtained for a constant minimum 
average volume, vl = 10 -29 m 3 , while the dashed curve is obtained for the constant maximum average 
volume, v’h = lCF 28 m 3 . The highest and lowest temperatures are given as 7 h = 177 K and T l = 16.6 K, 
respectively. 

A Q = jNkD(v h 2/3 - Vl 2/3 ) <0 (31) 

and 

(n 2 /2D)Tl(v 2 ^ - v 2/3 ) 

s =- 2 -7- L --. (32) 

r H ln(v H /v L ) - (3D/5 )(v“ 2/3 - v H 2/3 ) 


From Eq. (32) we can plot the ratio of the COPs versus the temperature ratio curves 
for three different volume ratios of two isochoric processes, as shown in Fig. 5, 
where the maximum average volume and the lowest temperature are given as v H = 
10™ 28 m 3 and Tf = 1 K, respectively. For a constant maximum average volume, 
vh = 10~ 28 m 3 , or a constant minimum average volume, vl = 10~ 29 m 3 , the ratio of 
the COPs versus the volume ratio curves are shown in Fig. 6, where the highest and 
lowest temperatures are given as 7 h = 177 K and T^—\ K, respectively. It is 
clearly seen from Fig. 6 that for a constant vl, the coefficient of performance 
increases with an increase of r v , while for a constant v H , the COP decreases with an 
increase of r v . These characteristics are similar to those in Fig. 4. 
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Fig. 5. The ratio of the COP of a Fermi Stirling refrigeration cycle to that of a reversible Carnot refrig¬ 
eration cycle versus the temperature ratio curves for three different volume ratios of two isochoric pro¬ 
cesses. The curves are presented for 3 He. The maximum average volume and lowest temperature are given 
as v’h = 10 _28 m 3 and 7 l = 1 K, respectively. 



Fig. 6. The ratio of the COP of a Fermi Stirling refrigeration cycle to that of a reversible Carnot refrig¬ 
eration cycle versus the volume ratio curves for 3 He. The solid curve is obtained for a constant minimum 
volume, vl = l(C 29 m 3 , while the dashed curve is obtained for a constant maximum average volume 
v H = 10~ 28 m 3 . The highest and lowest temperatures are given as 7 h = 177 K and If = 1 K, respectively. 
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5. Conclusions 

A Stirling refrigeration cycle working with an ideal Fermi gas does not possess the 
condition of perfect regeneration due to the quantum degeneracy of the Fermi gas. 
The COP of the cycle is smaller than that of a classical Stirling refrigeration cycle, 
whose COP is equal to that of a reversible Carnot refrigerator through the use of a 
reversible regenerator. Under the condition of strong gas-degeneracy, the COP of a 
Fermi Stirling refrigeration cycle, to a first approximation only depends on the 
temperatures of the heat reservoirs. Under the condition of weak gas degeneracy, 
the COP of the cycle is not only a function of temperatures but also dependent on 
the volumes and other parameters. 
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Appendix A 

In low temperature and high density conditions, i.e. under the condition of strong 
gas-degeneracy, the Fermi function f„(z) can be expressed as asymptotic expansions 
in powers of the quantity (Inc) -1 [13], i.e. 


/572(c) = 


15V^ 


(Inc) 


5/2 


1 +^-0nz) 2 


+ • 


(Al) 


fi/2(z) = w~7= (Inc) 

5^/71 


3/2 


TV 


l+l-(lnr)- 2 + 


(A2) 


Using Eqs. (Al ) and (A2), we can obtain the expression of the correction factor as 


2 1 + o n (l n “) 2 d- 

F(z) =f5/2{z)/fy2(z) = ~ ln(r)- =- 

1+-7r 2 (lnc) _2 H- 

O 


§ln(.-) 


Jt~ 


l+ T (lnc)- 2 +. 


(A3) 


To a first approximation [13], 
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1 £ f 

lm = tf 



(A4) 


( 3„ \ 2/3 (Ir \ 

where /x F = I-I I — is the Fermi energy. Substituting Eq. (A4) into (A3), to 

\47TgJ \2mJ 


a first approximation we can obtain 

2T f jt 2 T 

F( T, v) =- 1 -, 

V 5 T 6 T f 

where 7>(v) = D/v 2/3 is the Fermi temperature. 


(AS) 


Appendix B 


In high temperature and low-density conditions, i.e. under the condition of weak 
gas degeneracy, the Fermi function f,(z) can be expressed as an asymptotic expan¬ 
sion in powers of the quantity ( 2 ) [13], 


2 2 2 3 

/5/2(2) = - — 95/2 + 35/2 ' 


(Bl) 


2 2 2 3 
f3/ 2 (z') = 2 — 23/2 33/2 


(B2) 


Using Eq. (Bl) and (B2), we can obtain the expression of the correction factor as 


m = 


f5/2(2) 
f3/2(2) 


1 + 


-p2 + • • • . 

4V2 


(B3) 


Solving Eqs. (2) and (B2), we can obtain 

y =f3/i(z) = z~ 2 2 / 2 3/2 + 2 3 / 3 3/ '-, 


(B4) 


where 

)3 

y = - (B5) 

gv 

and 

2 = b\y + b 2 y 2 H-. (B 6 ) 

Substituting Eq. (B 6 ) into Eq. (B4) and comparing the coefficients of the same 
powers of y, we can obtain b\ — 1, b 2 — 1/2 3/2 ,.. Consequently, one can obtain 
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,_^ 3 1 ^ 3 V 

gv 2 3 / 2 Vgvy 


(B7) 


Substituting Eq. (B7) into Eq. (B3), to a first approximation we can obtain the 
expression for the correction factor as 

F(T, v) = 1 + E/(T y2 v). (B8) 

In the first approximation, Eq. (B7) can be written as 

z(T, v) = 4V2E/(T 3/2 v)[l + 2E/(T 3/2 v)]. (B9) 
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